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We discuss and calculate parity conserving (PC) and parity violating (PV) geo- 
metric phases for the metastable 2S states of hydrogen and deuterium. The atoms 
are supposed to be subjected to slowly varying electric and magnetic fields which act 
as external parameters for the atoms. Geometric flux density fields are introduced 
which allow for an easy overview how to choose the paths in parameter space in order 
to obtain only PC or only PV geometric phases. The PV phases are calculated in the 
Standard Model of particle physics. Even if numerically they come out small they 
have interest of principle as a new manifestation of parity violation in atomic physics. 
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1. Introduction 

In this article we discuss the atomic states of hydrogen and deuterium with principal 
quantum number n = 2 in slowly varying electric and magnetic fields. We are interested in 
parity (P) violating effects for the metastable 2S states. 

In the accompanying paper [1] we have already discussed the motivation for our study. 
There we have made a detailed investigation of the adiabatic limit for a system containing 

n 

metastable and short lived states. In the following the sections of [1| shall be quoted as 
section 1.1 etc. and equation numbers as (1. 1.1) etc. 

Our present paper II is organised as follows. In section [2] we review briefly the parity 
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violating Hamiltonian relevant for our atomic systems. We define our notation for the n = 2 
states and give the mass matrix for these states in external electric and magnetic fields. In 
this paper we always consider states at rest. The case of atoms travelling in the atomic 
beam interferometer [2] will be treated elsewhere. In section [3] we study the adiabatic limit 
for the metastable 2S states using the results of section Ij6] We identify the P-conserving 
(PC) and P- violating (PV) contributions to the geometric phases. In section H] we calculate 
these phases for the case that the geometric phase factor matrix is diagonal. In this case we 
have only abelian geometric phases. We give a graphical representation of these geometric 
phases in terms of surface integrals over geometric flux densities in the space of electric and 
magnetic fields. This allows us to see in an easy way how the paths in parameter space 
have to be chosen in order to get either pure PC or pure PV phases. Section [5] contains our 
conclusions and an outlook. In appendix |A] we collect tables giving the relevant numerical 
quantities of our systems. Appendix [B] contains a discussion of the adiabaticity condition 
for our concrete cases and appendix [0 gives the detailed calculations for the geometric flux 
densities. 

We use units h = c = 1 if other units are not explicitly indicated. 



2. Hamiltonian and state vectors 
2.1. The P-violating Hamiltonian 

In the framework of the Standard Model (SM) the effective P-violating Hamil- 

tonian relevant for atomic physics is due to the exchange of a Z boson between the atomic 
electrons and the quarks in the nucleus. In terms of the electron and quark field operators 
we have 

flpv = + (2-1) 
H<» = -j=f d Z x 2g\e{x)^ lh e{x) ^J2g q v q(x) llM q(x)j , (2.2) 

= ~ J d 3 x 2g e v e(x)^e(x) \^2g\q{x) llllb q{x) J , (2.3) 

where q = u,d,s, neglecting possible contributions from the heavy quarks c,b,t. The g e AV 
and g\ v are the neutral current coupling constants for the electron and the quarks, respec- 



tively, and G is Fermi's constant. For the notation and the definitions see 6] and [7]. In the 
framework of the SM the coupling constants of the weak neutral current are 



g e v = -\ + 2sin 2 ?V , 
9v = I - | sin2 $w , 
—\ + |sin 2 ?V , 
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(2.4) 
(2.5) 
(2.6) 



where dy/ is the weak mixing angle. 

For our study of atomic parity violation (APV) in light atoms, it is sufficient to consider 
a point-like, infinitely heavy, nucleus along with a nonrelativistic approximation of Hpy (for 
details see [7J]). In the nonrelativistic reduction, the effective P- violating Hamiltonians (12.21) . 
(12. 3p for an atom with proton number Z and neutron number N read 

G 1 



(i) 

PV 



Q$(Z,N) {5 3 (x)(*-p) + (<T-p)6 3 (x)} 



(2) 
PV 



G 1 



Q$(Z, N) {5 3 (x)(I ■ tr)(tr ■ p) + {tr ■ p)(I ■ <r)6 3 (x)} 



(2.7) 
(2.8) 



4\/2m e 

Here er and p are the Pauli spin matrix vector and the momentum operator for the electron, 

fl 2") 

respectively, and I is the nuclear spin operator. The Qy} (Z,N) are the weak charges of 
the atomic nucleus, given by 



Q$(Z, N) = -Ag\ {g"(2Z + N) + g d v {Z + 2N)} 
Q$(Z,N)= 4g e v J2g q Aj^(Z,N) . 



(2.9) 
(2.10) 



The weak charge Q^)(Z, N) exists, of course, only for nuclei with spin I ^ 0. The quantities 
Aq(Z, N) are the total polarisations of the nucleus carried by the quark species q. In the 
SM we have with to (1231) 



Q ( w\z,N) 



'l-4sm 2 $ w )Z - N 



(2.11) 

Q$(Z,N) = -4sin 2 ^) [Au(Z,N) - Ad(Z,N) - As{Z,N)} . (2.12) 

It was pointed out recently in [sj that the measurement of Q^(Z,N) is still the main 
motivation for APV studies, since it allows for a determination of sin 2 $yK, which is comple- 
mentary to high-energy physics. Such low-energy determinations of sin 2 i?w have, of course, 
been done with neutrino and electron scattering experiments. For a review see section 10 in 



19] 



9|. It is a challenge for experimentalists to see which method can give the best precision. As 
is clear from (12.121) APV measurements can also contribute to the study of the spin structure 
of light nuclei. For recent reviews of the available information on the spin structure of the 



nucleons see for instance 



2.2. The n = 2 states for hydrogen and deuterium 

We consider the subspace of atomic states with principal quantum number n = 2 for 
hydrogen and deuterium. As discussed in section IJ2], the effective Schrodinger equation for 
these systems reads in the Wigner-Weisskopf approximation 

i^\t) = JL{t)\t) , (2-13) 

where is the non-hermitian mass matrix and \t) the state vector of the undecayed atom 

in the (n = 2)-subspace. Here and in the following we use, in contrast to the companion 
paper I, the time t instead of the reduced time r (I J2.3H as variable for time dependent 
quantities. That is, we set tq = T in (I J2.3H . As in paper I we suppose that ^rff(t) is 
diagonalisable for all times t. Then we have for each time t a complete set of right and left 
eigenvectors satisfying 

Jt_{t)\a,t) = E(a,t)\a,t) , 

(a,t\jjf{t) = (a,t\E(a,t) , 

i (2-14) 
E{a,t) = E R {a,t) - -T{a,t) , 

(a = l,...,N) . 

Here E(a,t) are the complex energy eigenvalues which can be separated into a real part 
Eji(a,t) and an imaginary part — hT(a, t) leading to an exponential decay of the atomic 
state with decay rate F(a,t) (see section Ij2j) . 

For hydrogen with principal quantum number n = 2 there are N = 16 basis states, 
whereas for deuterium we have N = 24 states. As a set of basis vectors we choose coupled 
states of nuclear spin |/, ^3), electron spin I^S^) and electron orbital angular momentum 
\n, L, L3). We denote these basis states by \nLj, F, F3), where J is the total angular momen- 
tum of the electron and F, F 3 are the quantum numbers for the total angular momentum of 
the atom. 



5 



The mass matrix jtf. (t) contains contributions from the external fields and the PV Hamil- 
tonians (12.71) . (12.81) . The PV Hamiltonians mix the 25* and 2P states of hydrogen-like atoms 
in the (n = 2)-subspace. In the Coulomb approximation for the wave functions we get the 
following matrix elements 

(2S 1/2 ,F',F^i\2P l/2 ,F,F 3 ) = -i5 1 (Z,N)L(Z,N)8 FI>F 5 ntF3 , (2.15) 
(2S 1/2 , F', F" 3 \H$\2P 1/2 , F, F 3 ) = -iS 2 {Z, N)L(Z, N) [F(F + 1) - 1(1 + 1) - |] 

X 5 F r F S F > F3 , (2.16) 



I-l|<F,F'<I + i, . 



and 



(2S 1/2 ,F',F^\H^\2P 3/2 ,F,F 3 ) = , 

(2S 1/2 ,F',F^\H^\2P 3/2 ,F,F 3 ) = , (2.17) 



(|/-±|<F'</+±; |/-§|<F</+§) 
In (I2.15p . (12.161) define the PV parameters 

V3G Q$(Z,N) 



5i(Z,N) = - 



64nV2 r%(Z)m e L(Z,N) 

L(Z,N) 



-2.68827(3) aO-^feeV, ^ 



(i = l,2) . 

Here tb(Z) = (Zam^ 1 = Z~ l rs(l) is the first Bohr radius and L(Z,N) = E(2Si/ 2 ) — 
E(2Pi/ 2 ) is the Lamb shift for a hydrogen-like atom with proton number Z and neutron 
number N. The numerical factor in (I2.18P was obtained using the values for G, rs(l) and 
m e from 

With the Lamb shifts listed in table [3] of appendix [X] we get for ordinary hydrogen 

5,(1,0) = -6.14477(6) -10- 12 gSJ(l,0) , (z = l,2) , (2.19) 

and for deuterium 

5,(1, 1) = -6.13671(6) • lO" 12 Q$(l, 1) , (i = 1, 2) . (2.20) 

The numerical values of Si <2 for hydrogen and deuterium can be found in table [3] of appendix 
lAl The mass matrix for these atoms in the (n = 2)-subspace including PV contributions 
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can now be written as 

Mt) =^o-D- £(t) - fi ■ B(t) + 5x^\ + b 2 J^\ , (2.21) 

where we define for i = 1 , 2 

5iMv = (QL'j,, F', F* 3 \H®\2Lj, F, F 3 >) . (2.22) 

The mass matrix for zero external fields, 

JU?>Xi S 2 ) + + 5 2 ^pl , (2.23) 

and the matrix representations D and /x of the electric and magnetic dipole operators are 
shown in tables EJ [6] and [7] of appendix [A], respectively. 

At the end of this section, we introduce some notation. The atomic states for zero external 
fields, that is the eigenstates of ^_ (Si, S 2 ), will be written either in the short-hand notation 
|ck, or, with explicit quantum numbers, in the form \2Lj, F, F 3 ,5i,5 2 ). For non-zero 

electric and magnetic fields the eigenstates of the full mass matrix ^_{t) will be written in 
the forms \a, t, 61, 5 2 ), \a,£(t),B(t),5i,8 2 ), and \2Lj, F, F 3 , £(t), B(t), Si, 62), depending on 
the context. 



3. Adiabatic limit and geometric phases 

In the absence of external fields the 25* states are metastable, the 2P states decay fast. 
The theoretical lifetimes from recent QED calculations are 

r s = Tg 1 = 8.2207- 1 s = 0.1216 s , (3.1) 
Tp = TJ, 1 = (6.2649 ■ lO 8 )^ 1 s = 1.596 ■ 1(T 9 s , (3.2) 

see {14I An external electric field mixes the 2S and 2P states. The lifetimes of these 
mixed states are shown for hydrogen in figured! similar results can be obtained for deuterium. 
For the numbering scheme see table H] of appendix |A] Note that the metastable states are 
labelled in a different way as in I, where the numbering was a = 1, . . . , M for a total of M 
metastable states. 

In the electric field the 25" states, that is, the atomic states with a G {9, 10, 11, 12} = 
for hydrogen and a G {13, 14, 15, 16, 17, 18} = for deuterium quickly become more and 
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FIG. 1: The lifetimes of the mixed 2Si/ 2 and 2Pi/ 2 states of hydrogen in an external electric field 
S = £e 3 . Figure (a) shows the 2S , 1 / 2 states with labels a £ {9, 10, 11, 12}, figure (b) shows the 2P 1 / 2 states 
labelled by a £ {13, 14, 15, 16}. The lifetimes of the 2P 3 / 2 states with labels a £ {1, . . . , 8} are not shown 
here. The variation of the 2P 3 / 2 lifetimes is about two orders of magnitude smaller than that of the 2Piy 2 
lifetimes. For the numbering scheme see table |4] of appendix [A] 
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more unstable with increasing \£\. Still, for 

\£\ < 250V/cm (3.3) 
we have for the lifetimes of these states 

r a >10- 8 s, (ae.J% D ), (3.4) 

and 

rjrp > 5 (3.5) 

for a G Jm ,D an d P an y °f the IP states. The ratio of lifetimes ( 13.51) should be sufficient 
in order to apply the results of I to the systems considered here. In the following we thus 
have to keep in mind the restriction (13.31) in the discussion of adiabatic limits and geometric 
phases. 

We suppose that our atom is subjected to slowly varying electric and magnetic fields. In 
I we have — for mathematical convenience — worked with the limit T — > oo, making the 
variation of the fields smaller and smaller in real time. In reality we shall identify T with 
the observation time being of the order of t$ = T^ 1 . We discuss the resulting adiabaticity 
condition in appendix [Bj where based on the results of I we find that typically the rate of 
change of the external fields must satisfy 



1 

max — 

te[o,T] S 

1 

max — 

te[o,T] B 



dS(t) 



dt 
dB(t) 





1 


< 






rs 




1 


< 






TS 



(3.6) 

9t --• < 37 > 
where £ = 477.3 V/cm and B = 45.65 mT, see (fRlOlh (jBlBl . flR20|) and (IR22D . In the 
following we always suppose this to hold. 

Now we can apply the results of section IJ6J We suppose that we start at time t = with 
metastable states only, 

MO) = I (3.8) 
I otherwise . 

The time evolution of the metastable states is then given by (1 16. II) . neglecting non-adiabatic 
corrections, 

Mt) = J2 U °fiV)xf- (3-9) 
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B 3 [mT] B'i [mT] 

(a) (b) 

FIG. 2: The Breit-Rabi diagrams for the metastable states of hydrogen (a) and deuterium (b). The 
energy values are relative to the centre of zero-field energies of the 2P 1 / 2 states. 

Here and in the following we use J m always understanding it as for hydrogen and J® for 
deuterium. In the present paper we shall only discuss the case where, by suitable magnetic 
fields, the eigenenergies of the metastable states are separated such that condition (1 13. 2p 
holds. We have then, according to I, the case that the components of the metastable states 
evolve independently and no mixing occurs. In other words, the geometric phase matrix is 
diagonal (see (I J6.9f) -(I J6.11[) ) and we have (with r = T, r = t) 

Mt) = exp [-i<p a (t) + h aa (t)] (3.10) 

for a G J m . Here (p a (t) is the dynamical phase (see (I J3.4j) ) 

(p a (t) = [ dt' E(a,t') (3.11) 
Jo 

and jaait) is the geometric phase (1 13.51) given here by 

ft — Q 

j aa (t,5 1 ,6 2 ) = J dt (a,V,6i,6 2 \i-Qjj\a,1f,6i,62) ■ (3.12) 

In detail we shall suppose that there is a constant magnetic field B in 3-direction and 
that only the electric field £{t) varies with time. The magnetic field alone will lead to the 
usual Breit-Rabi diagram of the energy levels. This is shown in figures [2^ and [2b for the 
hydrogen and deuterium 2S states, respectively. 
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For a magnetic field \S\ > 1 mT the splittings in energy AE of the Zeeman levels satisfy 

\AE\ > AE = 0.02 /ieV , 

(3.13) 

AE /h = 5 MHz . 

Then the frequencies associated with the variation of the electric field should be much smaller 
than the one of (I3.13p . That is, we consider the Fourier analysis of the electric field strength: 

/OO J 
-e-^£(u). (3.14) 
■oo 27T 

The requirement is that ^ only for 

\u\ < \AE /H\ = 3- 10V 1 . (3.15) 

In other words, the variation of £(t) should be negligible over time intervals of order 
h/AEg = 2 • 10 _7 s. Then we will have a geometric phase matrix which is diagonal. 

In |l6j it was shown that time reversal invariance forbids shifts linear in the P-violating 
parameters for the eigenenergies of states in a spatially constant electromagnetic field. This 
implies, in particular, that the energies E(a,t) and the dynamical phases (p a {t) have no 
5-linear terms. They are insensitive to P- violation if we neglect terms of order 5 2 . The 
geometric phases, on the other hand, do have ^-linear terms as we shall see explicitly below. 

In the following it is convenient to discuss the geometric phases using the parameter space 
1Z, which in our case is the six dimensional space of electric and magnetic field strengths. 
The path C C 1Z along which the atom is led in 1Z is given by a time dependent parameter 
vector 

R(t) = (£ 1 (t),£ 2 (t),£ 3 (t),B 1 (t),B 2 (t),B 3 (t)) . (3.16) 

The time dependence of the atomic basis states (12.141) is only through R. That is, we can 

set 

| a, t, 5 U 8 2 ) = | a, R(t), 5 l , S 2 ) (3-17) 



and thus 



d_ 

ot 



a,t,5 1 ,5 2 ) = V R \a,R,5 1 ,5 2 ) . (3.18) 

R(t) at 
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Hence, the integral (13. 12[) can be transformed into a line-integral in parameter space, 

j aa (C,S 1 ,6 2 ) = dR ■ (a, R, 5x, S 2 \iV R \a, R, 61, 5 2 ) . (3.19) 



In section H] we shall consider closed paths C and transform ( 13.191) into a surface-integral. 
From this we will derive the geometric flux densities that will turn out to be extremely useful 
for visualisation purposes of geometric phases. 

However, numerical calculations of geometric phases for a given parametrised field confi- 
guration can be performed quite easily by using (I3.12p . Introducing the local matrix elements 
of the time derivative, 

^(R(t),5 x ,5 2 ) = (a,R{t)TSi,S 2 \^\(3,R(t),Si,S 2 ) , (3.20) 

the geometric phase (I3.12p reads 



Jaa(t,5 1 ,5 2 ) = i [ dt' ^ a (R(t'),Si,S 2 ) • (3-21) 

Jo 



(In the notation of I we have for r = T a aa = i3[ aa , see (I J4.3P ). 

In order to study PC and PV contributions to the geometric phases separately, we use 
perturbation theory to expand the atomic states in powers of the PV parameters 5i )2 . This 
is discussed in appendix We find with ( lC.36j) ff. for the matrix elements ( 13.201) for (3 = a 

@ aa (R(t), 5 U 5 2 ) = @ PC , aa m)) + ^ P \ aa (R(t)) + S 2 ^ aa (R(t)) + 0(5 2 ) , (3.22) 

where 0(S 2 ) is the short-hand notation for 0(S 2 , 6%, 5\5 2 ) and 

2pc**W)) = (a^(t)\^°\R(t)) , (3.23) 

(R(t)U PC (R(t)) 



- pv ' Qal 1 )] £^ \ E(a(°\ R(t)) - £7( 7 W, R(t)) 



+ 



(3.24) 



E(a(°\R(t))-E(j(°),R(t)) 
M\r, aj Wt)) = ( a ^\R(t)\^ Y \^ \R(t)) , (3.25) 
(* = 1,2) • 

Here, \a^°\ R(t)) and E(a<®, R(t)) are the eigenstates and eigenvalues of the mass matrix, 
obtained from (12.211) by setting Si = 5 2 = 0, 

JL{R{t),0) = ^_(R(t), Si, S 2 )\ Si=52=Q — — D ■ £(t) — /£ • B{t) . (3.26) 
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From (I3.2ip - (j3.25p we obtain the PC and PV contributions for the geometric phases, 

laa(t, 6 U S 2 ) = lP C , aa (t) + 8nFV,aa(t) + S 2j^, aa {t) + 0(5 2 ) , (3.27) 
lPC,aa{t) = 1 / dt' ^pc,aaW)) , (3-28) 

Jo 

l%, aa {t) = i fdt' ^ >aa (R(t')) , (* = 1,2) ■ (3.29) 
Jo 

4. Geometric flux densities 

In the following we always consider a closed path C in parameter space. 

The direct evaluation of fl3.27l) - fl3.29l) is the easiest way to obtain numerical values of 
geometric phases for a given path. The problem is to find suitable paths in parameter space 
that maximise PV or PC geometric phases while keeping the decay of metastable states 
to a minimum. In this section we will introduce geometric flux densities in order to get 
a visualisation and a better understanding of geometric phases. We will also discuss our 
numerical results for PV geometric phases at the end of the section. 

By using the generalised Stokes' theorem the integral (I3.19p . for a closed path C, can 
be transformed into a surface-integral in parameter space. The detailed calculation can be 
found in appendix ICl and gives (see (1C.30I) ) 

laa{C) = I j£(R) ■ df (£) + [ JI*J(R) ■ dfW + [ l£*){R) . (4.1) 

Here, T is a two-dimensional surface in the parameter space 1Z with C = dT . 

Thus, the geometric phases can be written as a sum of three integrals. The first inte- 
gral is evaluated in the space of electric field strengths, the second integral in the space 
of magnetic field strengths and the third in the full six-dimensional parameter space. The 
integrands J'^li.R) an d >?^{R) are three-dimensional vector fields in the spaces of elec- 
tric and magnetic field strengths, respectively. They can be interpreted as geometric flux 
densities and are a useful tool for visualisation and understanding of geometric phases as we 
will demonstrate below. 

The geometric flux densities can be split into PC and PV parts by using perturbation 
theory. The details of this calculation can be found in appendix IC.31 where we have used 
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only one combined PV parameter S (see ( 1C.31j) ). 

5= (Sl + 6 2 2 ) 1/2 . (4.2) 

From the PV geometric flux density in electric-field space (1C.48I) we can easily derive the 
nuclear spin independent and dependent flux densities. We obtain from flC.46l) - flC.48P the 
geometric flux densities in electric field space as 

Ji%(R, Si, h) = J%r\R) + 8xJ8™(R) + B*J$™(R) , (4.3) 

J, aa (*)-- 2 ^[2_, (E{ai0)jR) _ mo)jR)r ) » (44) 

r (£,pv) / jy\ _ i V- ( Dl^(R)Dilp a (R) + D™ ap (R)D™ a (R) - (i <-> j) 



V (E(a(°),R)-E(l3(°),R)y 



(4.5) 



for £ = 1,2, 3. Here the matrix elements of the electric dipole operator enter. For D p (R) 
see (IC.44p . The two PV contributions to the electric dipole operator follow from f]0.45|) and 
read 

^ ^ E(a(°\R) - E(<yW, R) ^ E(p<®, R) - S( 7 (°), i?) 1 J 

for x= 1,2. 

Similar results follow for the geometric flux densities in magnetic field space by making 
the replacements S — > B and .D — > /j, in (I4.3p - fl4.6l) . 

We will restrict ourselves in the following to the case of constant magnetic field. We 
are then dealing only with the three dimensional electric-field-strength parameter space. 
Correspondingly, the closed path C and the surface T with C = dT refer to this space. The 
geometric phase fl4.ll) reads then 



laa{C) = J^{R) ■ (for B = const) . (4.7) 

We will now study Sf^l{R) m some detail for the metastable states of hydrogen and 
deuterium. The numerical calculation of these vector fields is done for magnetic fields 
pointing in 3-direction. We show first some examples of flux density fields in electric-field 
space. In figures [3] and 0] the PC flux densities of the hydrogen states with a G {9, 10} 
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FIG. 3: The PC geometric flux density field of the hydrogen state |9) = \2Si/ 2} 1, 1, £, B = 1 mT) in electric-field space. The diagram on the l.h.s. 
shows the real parts of the flux density vectors, the diagram in the middle shows the imaginary parts. The diagram on the r.h.s. indicates the position 
of the plotted plane in electric field space. 




FIG. 4: The PC geometric flux density field of the hydrogen state 1 10) = |25i/ 2 , 1, 0,£ , B = 1 mT) in electric-field space. The diagram on the l.h.s. 
shows the real parts of the flux density vectors, the diagram in the middle shows the imaginary parts. The diagram on the r.h.s. indicates the position 
of the plotted plane in electric field space. 
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are shown. Rotational invariance of the fields around the ^-axis is indicated here, as the 
plots of the £i-£ 3 -plane are identical to the plots of the £ 2 -£ 3 -plane. In fig- E] a toroidal 
structure of the flux density field shows up, whereas in fig. H] a source and a sink can be 
seen. The corresponding plots for the other metastable hydrogen and deuterium states look 
very similar and can be found online fl7 |. 

Figure [5] shows the nuclear spin dependent PV flux density of the hydrogen state |9) = 
125*1/2, 1, 1, £, B) in the £i-£ 2 -plane fc> r £3 = 0. The solenoidal character and the cylindrical 
symmetry of the flux density can be seen. Many similar plots for different values of £3 and 
all other metastable hydrogen and deuterium states can also be found online 

From a study of many vector plots we find that {£, B) has only radial and 3- 

components whereas J~^aa\^^) ( x = 1>2) has only an azimuthal component. At least, 
this holds to the accuracy of the numerical calculations in all cases studied so far. 

Furthermore, the vector plots directly confirm the PC and PV character of the vector 
fields J^f c) (£,B) and J { ^\S,B), (x = 1,2). This is discussed in appendix [03 
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(IC.64I) and (1C.67I) . respectively. The plots also show that in all cases studied so far the PV 
vector fields in electric field space are source free, that is we have, to our numerical accuracy, 

V £ -JI% Y \S,B) = (4.8) 

for x = 1,2, whereas the PC flux fields can have sources and sinks. 

These vector field plots show that the PC and PV flux densities are orthogonal to each 
other and that they, furthermore, make it easy to choose closed paths in parameter space 
which maximise PC or PV geometric phases, respectively. Two examples of closed paths in 
parameter space that will give purely PV geometric phases are shown in fig. [6j The path 
(a) — consisting of two-point symmetric loops — gives a purely PV geometric phase for an 
arbitrary constant magnetic field since the PC contributions of the two loops cancel each 
other due to ( 1C.64I) . whereas due to ( 1C.67I) the PV contributions will be the same for each 
loop and add up. 

The PV geometric phases of all metastable hydrogen and deuterium states have been 
calculated numerically by evaluating the surface-integral (14.71) over the PV geometric flux 
densities (14. 5 p for a rectangular surface as shown in fig. [6b. In detail, we define the surface 

jr = {£. g x = 0) £ 2 e [0,250 V/cm] , £ 3 e [-125 V/cm, 125 V/cm]} (4.9) 
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FIG. 5: The nuclear spin dependent PV flux density field of the hydrogen state |9) = \2Si/ 2 , 1,1,£,B = 1 mT) in the £i-£2-plane with £ 3 = 0. The 
figure shows from left to right the real and imaginary part of the flux density and the location of the plane in electric field space. 
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(a) 



(b) 



So 



FIG. 6: Two example paths for a purely PV geometric phase, assuming a constant magnetic field in 
3-direction. 



and use a constant magnetic field 



B = (0,0, lmT) . 



(4.10) 



We have also evaluated the line-integrals (13.291) along dT numerically in order to cross-check 
our results which are listed in table [TJ We find that the results from both methods match 
nicely and give an indication of the accuracy of the numerical calculations. The total PV 
geometric phases are (see (13.271) ) 



Tpv.i 



(4.11) 



with the PV parameters 5\ t 2 taken from table [3j The 7pv,»a are listed in table [2j 

We will now discuss the results shown in table |2j First of all, the metastable hydrogen 
states |9) and 1 10) have the largest PV geometric phases. In figured the corresponding real 
parts of the nuclear spin dependent fluxes tfiaa ' df^/\df^\ are shown. They 

represent the dominant contributions to 7pv,oa for these states. The areas of the largest 
contributions to the flux are centred around £3 = and extend to higher values of £2- In 
fact, the PV geometric flux can be increased by one order of magnitude for a surface with 
£2 £ [0, 1000 V/cm], but then, of course, (I3.3P no longer holds for every point on that surface 
and the mixed 2S states can in general no longer be considered as metastable, see appendix 

m 

The quoted uncertainties in the results of table [2] only take into account the errors of 
the PV parameters di^ given in table [3j The numerical uncertainties have been neglected 
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a 


\2S 1/2 ,F,F 3 ) 


7p\ ; 

surface 


r ,aa 

line 


Tp\ 

surface 


r ,aa 

line 


9 

10 
11 


|2S\/ 2 ,M) 
|2^i/2,l,0) 
|25 1/2 ,1,-1) 


-0.1966 + 0.0290i 
0.1438 - 0.0194i 
0.0394 - 0.0078i 


-0.1918 + 0.0287i 
0.1438 - 0.0194i 
0.0352 - 0.0074i 


-3.0547 + 0.3101i 
4.2083 - 0.3872i 
-0.9889 + 0.0584i 


-3.0548 + 0.3099i 
4.2132 - 0.3873i 
-0.9936 + 0.0586i 


12 


|2^i/2,0,0) 


0.0339 - 0.0055i 


0.0330 - 0.0054i 


-0.1610 + 0.0179i 


-0.1610 + 0.0180i 



(a) 



a 


\2S 1/2 ,F,F 3 ) 


surface 


',oa 

line 


(2; 
Tp< 

surface 


r ,aa 

line 


13 
14 
15 
16 


|2<^i/ 2 , §, |) 
25*1/2, f , 2) 

l 2 ^!^, §, -5) 

25*1/2, §, -|) 


-0.1382 + 0.0218i 
-0.0767 + 0.0119i 
-0.0123 + 0.0019i 
0.0924 - 0.0149i 


-0.1333 + 0.0214i 
-0.0733 + 0.0116i 
-0.0108 + 0.0018i 
0.0876 - 0.0146i 


-20.004 + 2.247i 
11.642 - 1.296i 
13.441 - 1.456i 
-0.235 - 0.040i 


-20.033 + 2.249i 
11.665 - 1.297i 
13.460 - 1.458i 
-0.242 - 0.040i 


17 
18 


25*1/2, 2J 2) 
25i/ 2 , 2", -f) 


0.1284 - 0.0200i 
0.0365 - 0.0059i 


0.1246 - 0.0197i 
0.0348 - 0.0058i 


-11.545 + 1.346i 
6.703 - 0.801i 


-11.564 + 1.348i 
6.715 - 0.802i 



(b) 



TABLE 1: Numerical results of PV geometric phases for metastable hydrogen (a) and deuterium (b). 

but can be estimated roughly by comparison of the results for the surface and the line- 
integration and seem to be around a few percent. The only PV geometric phase for which 
the surface- and the line-integration results do not coincide within la is the result for the 
metastable deuterium state 1 16) = |2S*i/ 2 , §, — §,£,25). The small uncertainty due to <5i j2 in 

7pv,i616 is due to its small nuclear spin dependent contribution 7 P vi 6 i 6 (see tabled]). For 

(2) 

all other hydrogen and deuterium states the contribution from 7 P y aa dominates over that 
from 7p VQ , Q in 7pv,aa- This, together with the fact that <5 2 for deuterium is consistent with 
zero and has a statistical error larger than 500% (see table [3]), is the reason for the large 
uncertainties due to 5i t2 (up to 50%) of all metastable deuterium states with a 7^ 16. 

In fig. [8] we show the real parts of the PV flux densities for the deuterium states 
1 13) = |2S*i/ 2 , §,§,£, 25) and 1 16) = \2Sy 2 , §,—§,£, 25). Figures [8k-c are similar to the 
corresponding hydrogen flux densities shown in fig. [71 Only fig. [HU, which shows the nu- 
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|2S\ /2 , F,F 3 ) 


7PV,ao [1 

surface 


0~ 13 rad] 

line 


9 

10 
11 


|2«S 1/2 ,1,1) 
|25 1/2 ,1,0) 

12^/2,1,-1) 


-20.96(59) +2.10(6)i 
29.23(81) - 2.67(7)i 
-7.07(19) +0.41(l)i 


-20.97(59) +2.10(6)i 
29.26(81) - 2.67(7)i 
-7.09(19) + 0.43(l)i 


12 


|2S\ /2 ,0,0) 


-1.23(3) + 0.141(3)i 


-1.23(3) +0.141(3)i 



(a) 











7PV,«a [10 13 rad] 
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2S1/2, 


F,F 3 ) 


surface 


line 


13 


25' 1 /2, 


3 

2" 


1) 


-8.7(3.4) + 1.34(38)i 


-8.4(3.4) + 1.32(38)i 


14 


25' 1 /2, 


3 
2- 


1) 


-4.1(2.0) +0.66(22)i 


-3.9(2.0) +0.64(22)i 


15 


|25i/2, 


3 

2' 


-1) 


-0.32(2.3) +0.07(24)i 


-0.23(2.3) +0.06(25)i 


16 


1 2^1/2 5 


3 

2' 


"I) 


5.404(36) - 0.876(8)i 


5.128(38) -0.856(7)i 


17 




1 

2' 


1) 


7.2(2.0) - 1.13(23)i 


7.0(2.0) - 1.12(23)i 


18 


25*1/2, 


1 

2" 


-1) 


2.3(1.1) -0.37(13)i 


2.2(1.1) — 0.36(14)1 



(b) 



TABLE 2: The total PV geometric phases f|4. 1 1[) for the metastable states of hydrogen (a) and deuterium 
(b). The errors quoted are discussed in the text. 

clear spin dependent PV flux density of |16), looks different. Compared to fig. [8b one can 
see a difference of two orders of magnitude between both flux densities. Furthermore, fig. 
[HU shows a change of sign between the maximum and minimum of the plotted flux density 
values. The geometric flux density also increases both in positive and negative ^-direction. 

5. Conclusions and Outlook 

In this work we have shown that PV geometric phases in the lightest metastable atoms, 
25* hydrogen and deuterium, exist. We have considered adiabatic changes of the electric field 
along suitable closed paths C with a constant magnetic field in order to avoid degeneracies 
of atomic energy levels. The metastable states of hydrogen and deuterium indeed pick up 
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FIG. 7: The real parts of the projected, nuclear spin dependent flux densities 

J { 2 £ a ^ V \£,B) ■ d/ (£) /|df (£) l of the metastable hydrogen states |9) (a) and 1 10) (b). Brighter areas 
correspond to a larger absolute value of the geometric flux. The solid black lines are contour lines for the 
levels shown on the right hand sides of the density plots. 



non-zero PV geometric phases in this way. 

The geometric flux densities provide us with a good way to visualise geometric phases as 
fluxes through areas T in parameter space with C = OF. In this way we could see how to 
choose the path C in order to get only a PC or only a PV geometric phase. 

We have then studied numerical values of geometric phases for a closed rectangular loop 
in electric-field space (see fig. Eb) where only PV and no PC phases occur. The PV phases 
have been calculated numerically both from the surface and from the line-integrals, see (j4.7p 
and (13.291) . respectively. The results agreed within the numerical uncertainty. The resulting 
PV geometric phases are listed in tables [I] and [2j The geometric contributions 7py aa can be 
large, see table [U The smallness of the results for the total PV geometric phase in table [2] 
is solely due to the smallness of the PV parameters 8\^- It remains to be seen if such small 
effects can be enhanced by suitable iterations or other means in order to make them large 
enough to be measurable. 

Even if the results for PV geometric phases presented in this paper are numerically small 
we find that they have an interest of principle. They show that geometric phases can carry 
information from parity violation in atoms. This was first discussed in 18] and has now 
been investigated in detail for the 2S states of hydrogen and deuterium. 

There are numerous things left to be studied in the future. So far, we have only inves- 
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FIG. 8: The real parts of the two PV contributions to the geometric fluxes of the metastable deuterium 
states 1 13) = |2S 1/2 , §, § , £, B) (figures (a), (b)) and 1 16) = |2Si/2, §, -§,£,B) (figures (c), (d)). In (a) 
and (c) the nuclear spin independent, whereas in (b) and (d) the nuclear spin dependent contributions are 
shown. 



tigated abelian geometric phases for constant magnetic fields. The investigation of abelian 
geometric phases for varying magnetic field strengths is a challenging task for further stud- 
ies. Furthermore, non-abelian geometric phases could reveal interesting new features of PV 
geometric phases. 
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Appendix A Values for quantities related to n = 2 hydrogen and deu- 
terium 

In this appendix we collect the numerical values for the quantities entering our calcula- 
tions for hydrogen and deuterium states with principal quantum number n — 2. We specify 
our numbering scheme for these states and give the expressions for the mass matrices at 
zero external fields, and for the electric and the magnetic dipole operators. 

In table [3] we present the numerical values for the weak charges Qyft , x = 1,2, the 
quantities Aq, the Lamb shift L = E(2Sy 2 ) — E(2Pi/ 2 ), and the fine structure splitting 
A = E(2P 3 / 2 ) — E(2Pi/2). The ground state hyperfine splitting energy is denoted by A. 
We have A = E(1S 1/2 ,F = 1) - E(1S 1/2 ,F = 0) for hydrogen and A = E(1S 1/2 ,F = 
3/2) — E(lSi/ 2 ,F = 1/2) for deuterium. The n = 2 states of hydrogen and deuterium 
in the absence of P- violation and for zero external field are denoted by \2Lj,F, F 3 ), where 
L, J, F and F 3 are the quantum numbers of the electron's orbital angular momentum, its 
total angular momentum, the total atomic angular momentum and its third component, 
respectively. The quantum numbers S for the electron spin and / for the nuclear spin are 
omitted, since these are fixed quantities for each atomic species. In the following the ordering 
of the atomic states in the matrix representations of operators is according to decreasing F 3 , 
F, J and L. 

In table H] we give the numbering scheme for the states which we consider. For electric 
field S and magnetic field 3 equal to zero we have the free 25' and 2P states. We write L, 
P, S since these states include the parity mixing due to Hpy (12.11) . 

Consider first atoms in a constant B-field pointing in positive 3-direction, 

B = Be 3 , B > . (A.l) 

The corresponding states \2Lj, F, F 3 , 0, Be 3 ) are obtained from those at B = by continu- 
ously turning on 3 in the form flA.ll) . Of course, for \3\ 7^ 0, F is no longer a good quantum 
number. Here it is merely a label for the states. 

We now choose a reference field S re f = B ve {e 3l £> re f > 0, below the first crossings in the 
Breit-Rabi diagrams, for instance £> re f = 0.05 mT. We define the states \2Lj, F, F 3 , £, 3) 
for arbitrary 3 fields in the neighbourhood of S re f as the states obtained continuously from 
\2Lj, F, F 3 , 0, re f) by turning on £ and (3 — B rc f) as X£ and X{3 — B rc f), respectively, with 
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{H 
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Ref. 


z 


1 


1 




N 





1 




I 


l 

2 


1 




Q$(Z,N) 


0.04532(64) 


-0.95468(64) 


fl2.ll} 


St{Z,N) 


-2.78(4) • 10~ 13 


58.59(4) • 10~ 13 


(|2.19|).(|2.20p 


Au(Z,N) - Ad(Z, N) 


1.2695(29) 





[9J 


As(Z,N) 


0.006(29) (7) 


0.012(58)(14) 


[19] 




-0.1145(31) 


0.0005(26) 


fl2.12j) 


HZ,N) 


7.04(19) • 10~ 13 


0.03(17) • 10~ 13 


(I2.l9|).(|2.20l) 


L(Z,N)/h 


1057.8440(24) MHz 


1059.2330(26) MHz 


[20] 


A(Z, N)/h 


10969.0416(48) MHz 


10972.0355(48) MHz 


[20J 


A(Z, N)/h 


1420.405751768(1) MHz 


327.384352522(2) MHz 


[2JJ 



TABLE 3: Values of parameters for numerical calculations. The weak mixing angle in the low energy 
limit, sin 2 = 0.23867(16), was taken from [22|. The uncertainty in 8\ is dominated by the uncertainty 

of sin 2 "&w The values of and 62 for deuterium are consistent with zero, according to the value of As. 

(2) 

The uncertainties in Q w and 62 for hydrogen are resulting from the errors of the weak mixing angle and 
As in equal shares. 

A G [0, 1]. Then both F and -F3 are only labels for these states and no longer good quantum 
numbers. 
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TABLE 4: The numbering scheme for the atomic states of hydrogen and deuterium. 
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Tables [5j M and [7J show the non-zero parts of the mass matrix for zero external fields, 
of the electric dipole operator .D and of the magnetic dipole operator \x for the n = 2 states 
of hydrogen. 

In tables M and [7] we use the spherical unit vectors, which are defined as 

e = e 3 , e ± = =f-^= (ex ± ie 2 ) , (A.2) 
where e« (z = 1,2,3) are the Cartesian unit vectors. For e±, the following relation holds: 

e% = -e T . (A.3) 



TABLE 5: The mass matrix n(<5i , S2) (|2.23| for the n = 2 states of hydrogen for the case of zero 
external fields. For the explanation of the variables A, L and A and their numerical values see the 
introduction of this appendix and table [3l The PV parameters 5\ t 2 can also be found in table [3l the decay 
rates Tp.s are given in (|3.ip and (|3.2p . 





2P3/2 


2,2 


2^3/2,2,1 


2P 3 /2,M 


2Pi/ 2 ,M 


25i/ 2 , 1, 1 


2^3/2,2,2 


A + T6T> 

-|r P 














2^3/2,2,1 







A + 4 

-|r P 











2P 3 /2,M 










A-^ 

. 96 


.4 

192^ 





2Pi/ 2 ,l,l 










.4 

192^ 


.4 i r 
96 - 2 L P 


i8 x L 


261/2, 1, 1 















T 4- ^ 
^ + 32 

-ITs 



(Table [5k) 





2P 3 /2,2,0 


2^3/2,1,0 


2Pi /2 ,l,0 


25 1/2 ,1,0 


2Pi/ 2 ,0,0 


2^3/2,0,0 


2^3/2,2,0 


^ y i6o 
-§r P 

















2^3/2,1,0 







A 











2Pi/ 2 ,l,0 





.4 

192^ 


96 ~ S 1 ^ 


iSxL 








2S 1/2 ,1,0 








-i<5 x L 


T -I- A 
L+ 32 








2Pi/ 2 ,0,0 














A i -p 
32 2 1 p 


+§i<5 2 L 


2S 1/2 ,0,0 














-iS x L 


j 3A 

L ~ 32 



(Table [5b) 





2^3/2,2,-1 


2P 3 /2,1,-1 


2^1/2,1,-1 


251/2, 1, -1 


2P 3 /2,2,-2 


2P3/2 


2, 


-1 


A + T6l) 














2P3/2 


1, 


-1 







A 
192^ 








2Pl/2 


1, 


-1 





A 
192^ 


A i T 
96 - 2 L P 


15 X L 
+^8 2 L 







1, 


-1 








-i<5iL 


^+32 

-Ps 





2^3/2 


2, 


-2 














A + 1M 



(Table [5b) 



TABLE 6: The suitably normalised electric dipole operator D/(e r»(111 for the n = 2 states of hydrogen. 





2^3/2, 2, 2 


2^3/2, 2, 1 


2P 3 /2,M 


2Pi/ 2 ,l,l 


25' 1 /2, 1, 1 


2^3/2,2,2 














-3e_ 


2^3/2,2,1 
















2P 3 /2,M 














-\/i e o 


2^1/2,1,1 














-\/3e 


25 1/2 ,1,1 


3e + 


3 

75 e ° 


-\fi e o 


-\/Se 






(Table W) 





2P 3 /2,2,0 


2^3/2,1,0 


2P 1/2 ,1,0 


251/2,1,0 


2Pi/ 2 ,0,0 


251/2,0,0 


2^3/2,2,1 



















2P 3 /2,M 











— »Ae 

V 2 





-^/6e_ 


2^1/2,1,1 


















2^/2, 1,1 


a/^ 

V 2 


V 2 


-\/3e_ 





\/3e_ 






(Table [6b) 



to 

CO 







2P 3 /2,2,1 


2P 3/2 ,1,1 


2Pi/ 2 ,l,l 


2^/2, 1,1 






2P 3/2 ,2,0 

















2P 3 / 2 ,1,0 











V 






2Pi/ 2 ,l,0 

















2^/2, 1,0 






V3e + 









2Pi/ 2 ,0,0 

















25 1/2 ,0,0 





\/6e + 


-V3e + 













(Table [6b) 










2P 3 /2,2,0 


2P 3/2 ,1,0 


2Pi/ 2 ,l,0 


2Si/ 2) l,0 


2Pi/ 2 ,0,0 


25i /2 ,0,0 


2^3/2,2,0 











\/6e 








2^3/2,1,0 

















\/6eo 


2Pi/ 2 ,l,0 

















-\/3e 


2S 1/2 ,1,0 


VQe 











-A/3e 





2Pi/ 2 ,0,0 











-\/3e 








2Si/ 2 ,0,0 





\/6e 


-V3e 












(Table EH) 





2^3/2,2,-1 


2^3/2,1,-1 


2Pi/ 2 ,l,-l 


251/2, 1, -1 


2^3/2,2,0 











V 2 


2ify a ,l,0 











V 2 e - 


2Pi/ 2 ,l,0 











-\/3e_ 


25 1/2 ,1,0 


V2 e ~ 


V 2 







2P 1/2 ,0,0 













25i /2 ,0,0 





^6e_ 


-\/3e_ 






(Table [6b) 





2P 3 /2,2,0 


2P 3 /2,1,0 


2Pi/ 2 ,l,0 


251/2,1,0 


2Pi/ 2 ,0,0 


251/2,0,0 


2-P3/2, 2, 


-1 



















2^3/2,1, 


-1 











\fi e + 





-y/Ee+ 


2^1/2,1, 


-1 











V3e + 







2-51/2, 1, 


-1 




7f e + 


V3e + 





V3e+ 






(Table [6f ) 





2P 3/2 ,2,-l 2P 3/2) 1,-1 2P 1/2 ,1,-1 2S 1/2 ,1,-1 


2P 3 /2,2,-2 


2 ^3/2 > 2 , -1 

2^3/2,1,-1 

2^1/2,1,-1 
2'5 1 /2, 1, -1 


^e 
OOO yjeo 
\/3e 
^e \[\ e o V3e 





3e 


2P 3 /2,2,-2 


-3e+ 






(Table Efe) 



TABLE 7: The suitably normalised magnetic dipole operator /jl/^b for the n = 2 states of hydrogen, where \xb = eh/(2m e ) is the Bohr magneton and 
g = 2.002319304(76) is the Lande factor of the electron 





2^3/2,2,2 


2^3/2,2,1 


2^3/2,1,1 


2^1/2,1,1 


25 V 2,1,1 


2^3/2 


2,2 




^(9+2) 
4 e 


^(9+2) 
4^ 







2^3/2 


2,1 


4 e + 











2^3/2 


1,1 


^(9+2) 


4^ 




§73 e ° 





2-Pl/2 


1,1 






3^ 







25*1/2, 


1,1 














-fee 



(Table [7» 





2^3/2,2,0 


2^3/2,1,0 


2^1/2,1,0 


25*1/2,1,0 


2Pi/2,0,0 


25i /2 ,0,0 


2P 3 /2,2,1 

2^3/2,1,1 

2Pi/ 2 ,l,l 


^(9+2) 
4 e 

9+2 _ 
12 

6 e 


9+2 
4^ 

5(9+2) 
12 e - 

V^(ff-1) _ 
6 e 


6 e 
6 e - 









V2(g-l) 
3 e 

6 e 







25i/ 2 , 1, 1 











-9-e 

2 e - 





9-e 



(Table [7b) 



CO 

to 





2P 3/2 ,2,1 2P 3/2 ,1,1 2P 1/2 ,1,1 25i /2 ,l,l 


2^3/2,2,0 

2^3/2,1,0 

2Pi/ 2 ,l,0 
2S 1/2 ,1,0 
2Pi/ 2 ,0,0 
25 1/2 ,0,0 


4 L + 12 + 6 + 

.9+2 5( 9 +2) V^(S-1)_ n 
4^ + 12 L + 6 e + U 

9-1 V^(9-l) 9-4 

fe + 
^fV f 6 4 e + 
-fe+ 



(Table Eb) 





2^3/2,2,0 


2^3/2,1,0 


2P 1/2 ,1,0 


2S 1/2 ,1,0 


2Pi/ 2 ,0,0 


2^1/2,0,0 


2^3/2 


2,0 





-^o 













2 -P 3 /2 


1,0 


















2Pi/ 2 


1,0 


















2Si/2 5 


1,0 

















-feo 


2Pi/ 2 


0,0 







9- 4 „ 











25'i/ 2 , 


0,0 











-f e o 









(Table EH) 



2^3/2,2, 



1 2P 3/2 ,1,-1 2P 1/2j l,-l 25 1/2) 1,-1 



1/2) 



'1/2; 




2^3/2,2,0 
2^3/2,1,0 

2Pi/ 2 ,l,0 
25 1/2 ,1,0 
2Pi/ 2 ,0,0 
25i /2 ,0,0 



(Table 03) 





2P 3 /2,2,0 


2P 3 /2,1,0 


2Pi/ 2 ,l,0 


2Si /2 ,l,0 


2Pi/ 2 ,0,0 


25 1/2 ,0,0 


2^3/2) 2, 


-1 


V3( S +2) 
4 e + 


4V3 + 













2P 3 /2,1, 


-1 


12 e + 


5(5+2) 
12 e + 


6 + 





3 e + 





2Pi/ 2 ,l, 


-1 


6 + 


V5(fl-i) r 

6 + 


9-4 „ 

-V e + 





-V e + 





2»S'i/2, 1, 


-1 











|e + 





2 + 



(Table m) 





2P 3/2 ,2,-l 2P3 /2 ,1,-1 2P 1/2 ,1,-1 2S 1/2 ,1,-1 


2P 3 /2,2,-2 


2 ^3/2 > 2 , -1 

2^3/2,1,-1 

2^1/2,1,-1 
2'5 1 /2, 1, -1 


^en -^en -^en 
4 e ° 4V3 £ ° V6 £ ° U 

-S5fe "§73*0 

"W e ° -f73 e o 
fe 


^(9+2) 
4 e 

^(9+2) 
4^ £ - 




2P 3 /2,2,-2 


^(9+2) ^(9+2) 3-1 n 
4 C ^ 4^ £ + + 


2+g „ 



(Table Wg) 



36 



Tables El [9] and [TO] show the mass matrix ^ for zero external fields, the electric dipole 
operator D and the magnetic dipole operator fi for the n = 2 states of deuterium. 

TABLE 8: The mass matrix ^0(61,62) (|2.23[) for the n = 2 states of deuterium for the case of zero 
external fields. For the explanation of the variables A, L and A and their numerical values see the 
introduction of this appendix and tabled The PV parameters 5\ t 2 can also be found in table 02 the decay 
rates Tp.s are given in (|3.ip and (|3.2p . 





2-P3/2 


5 5 

2' 2 


2-P3/2 > 2 ' 1 


op 3 3 
Z/ 3/2) 2> 2 


op 3 3 
-^1/2) 2' 2 


00 3 3 
^1/2) 2' 2 


2-^3/2) 


5 

2' 


5 
2 


^ T 120 














2-P3/2 ! 


5 

2' 


3 
2 







-|r P 











2-P3/2) 


3 

2" 


3 
2 










" 180 


576 





2^1/2, 


3 

2' 


3 
2 










y/EA 
576 


72 ~ 2 L P 


i8 x L 
+i8 2 L 




3 

2' 


3 
2 













-i<5iL 
-i(5 2 I/ 





(Table [Si) 





2-P3/2> |; 2 


2^3/2, f > I 


2-Pl/2, f > \ 


25*1/2 1 f > | 


op 11 

^3/2 > 2> 2 


op 11 

*-TL/2) 2' 2 


oc 11 
z,5 l/2> 2' 2 


9 P , 
ZJ 3/2 


5 1 
2' 2 


A + — 

^ y 120 

-7,r P 




















^ r 3/2i 


3 1 
2' 2 





A A 

A — — 

180 


576 














2Pi/ 2 


3 1 
2' 2 





576 


A i -n 
72 _ 2 L P 


1U1 -LJ 

+i5 2 L 











2-S1/2 


3 1 
2' 2 








-iSxL 
-\b 2 L 


L + A 











2^3/2 


1 1 

2' 2 














A-A 

-ir P 


A 

288V2 





2 -Pl/2 


1 1 

2' 2 














.a 

288^ 


A 

36 

-|r P 


S6 X L 
-2i5 2 L 


25*1/2 


1 1 

2' 2 

















- vSxL 
+2\5 2 L 


T A 

h 12 



(Table [8jo) 





2-P3/2, f, _ 2 


2A/2, |, -3 


2^1/2, 3, -3 


25i/ 2 , §, -3 


op 1 1 

•^3/2, 3> -3 


op 1 1 

^1/2; 3> -3 


25i/ 2 , 3, -5 


9 P , 
■"3/2- 


5 

2' 


1 

2 


^ y 120 






















3 
2' 


1 

2 





A A 

-3^P 


VEA 
576 














2Pi/ 2 


3 
2' 


1 

2 





576 


72 - 2 L P 


1U1 -L^ 











25*1/2, 


3 
2- 


1 

2 








-iSxL 
— i5 2 L 


^ + 4 











2P 3 /2 


1 

2' 


1 

2 














-3^P 


.4 

288x72 





2^1/2 


1 

2' 


1 

2 














A 
288,72 


A 

36 

-|r P 


-2\5 2 L 


25i/ 2 , 


1 

2> 


1 

2 

















- WiZ 
+2i5 2 L 


T A 

L ~ T3 



(Table lb) 





2^3/2, 2' ~l 


OP 3 3 
^3/2 > 2 ' ~~ 2 


2 -Pl/2) |> ~i 


c 3 3 
^1/2) 2' ~2 


2-P3/2, 2' _ 2 


2^3/2, 


5 

2' 


3 
2 


-|r P 














2P3/2, 


3 
2' 


3 
2 





A — ^_ 

. 180 


576 








2P1/2, 


3 

2' 


3 
2 





576 


.4 i -n 
72 ~ 2 L P 


iSxL 
+i6 2 L 





25' 1 /2, 


3 

2' 


3 
2 








—id 2 L 


-frs 





2^3/2, 


5 

2' 


5 
2 














^ \ 120 



(Table EH) 



TABLE 9: The suitably normalised electric dipole operator D /(errt(l)) for the n — 2 states of deuterium. 





op 5 5 
-"3/2 j 2' 2 


2-P3/2 5 3' 1 


2-P3/2, |> | 


2pL/2> f , 1 


c 3 3 
^1/2 > 2» 2 


2P3/2, 


5 

2' 


5 
2 














-3e_ 


2P3/2, 


5 

2 ' 


3 
2 
















2P3/2, 


3 
2' 


3 
2 
















2Pl/2, 


3 

2 - 


3 
2 














— v/3e 


25' 1 /2, 


3 
2' 


3 
2 


3e + 




-2^0 


-\/3e 






(Table Eh) 



2-P3/2)|!3 2P 3 / 2 ,§, 3 2P 1 / 2 ,|,3 251/2, |, 3 2P 3 / 2 ,2,2" 2Pi/ 2 ,5)3 251/2, 3, 3 



3 1 



3 1 



3 1 



1 1 



1 1 



1 1 




(Table ©)) 



op 5 3 op , 3 3 op 3 3 no , 3 3 
Z/ 3/2)2'2 ^3/2) 2' 2 Z/ l/2> 2' 2 ^1/2) 2' 2 




(Table [9b) 




(Table [Si) 




(Table [9b) 



2-F3/2) 2' 2 



2^3/2, f , 2 



2^1/2, f , \ 



25i/ 2 , |, ^ 



2P, 



"3/2, I' I 



2P, 



1/2, 2> I 



251/2, 2' 2" 



2^3/2, |, 
2^3/2, |, 
2 -Pi/2, §, 
251/2, |, 



2P. 



3/2, 2' 
2^1/2, 2", 

25'i/2, ^, 







_3_ 
10' 

























§ e+ 





7§ e + 







A4 
V 3 e + 





































\/i e + 







2\/|e + 








(Table W) 



2-^3/2) 2P 3 / 2 ,|, 



2 2P 1/2 ,f 




(Table Efe) 





2-P3/2, I' ~l 


2-P3/2, f , -f 


2 -Pl/2, f , -| 


n c 3 3 
^1/2) 2' ~2 


2P3/2, 


5 

2' 


1 
2 











^ e 


2P3/2, 


3 
2' 


1 
2 











V 5 e ~ 


2^1/2 ) 


3 
2' 


1 
2 











-^2e_ 


2 <Si/2> 


3 

2' 


1 
2 




V 5 e ~ 


-^2e_ 





2P 3 /2, 


1 

2' 


1 
2 











V2 B - 


2Pi/ 2 , 


1 

2' 


1 
2 











-2e_ 


25'i/2) 


1 

2' 


1 
2 





^/5e_ 


-2e_ 






(Table [9fr) 





2-^3/2) §i - I 2 ^3/2) |> _ 2 


2 A/2> |) ~2 


25*1/2, §,- 5 2-^3/2, |> - 2 


2 A/2, ^> ~2 


25i/ 2 , 2> _ 2 


2^3/2, 


5 

2' 


3 
2 








- 3 7f e + 








2P3/2, 


3 
2' 


3 
2 








/f £ + 





— v/5e + 


2^1/2- 


3 
2' 


3 
2 








^/2e+ 





2e+ 


25*1/2, 


3 
2- 


3 
2 


io e + \/t e + 


^2e+ 


75 e + 


2e + 






(Table IB) 





2-P3/2, 25 ~| 


2^3/2, | ; -§ 


2-Pl/2> |> -§ 


c 3 3 
^1/2) 2' ~2 


op 5 5 
<^3/2) 2' ~2 


2P3/2, 


5 

2' 


3 
2 
















2^3/2, 


3 

2' 


3 
2 
















2 -Pl/2; 


3 

2' 


3 
2 
















25i/ 2 , 


3 

2' 


3 
2 


S^/feo 




\/3e 





3e+ 


2^3/2, 


5 

2' 


5 
2 











-3e_ 






(Table E|) 



TABLE 10: The suitably normalised magnetic dipole operator /J.//J.B for the n = 2 states of deuterium, where hb — eh/(2m e ) is the Bohr magneton 
and g = 2.002319304(76) is the Lande factor of the electron [23 1. 





op 5 5 
^3/2 > 2 ' 2 


od 53 oD 33 od 33 oc 33 
^3/2, 5, 5 ^3/2 > 2' 2 2 - r l/2>2'2 Z,5 l/2> 2' 2 


2^3/2 > !> | 


2+9 „ 
2^ e 


g+2 c g+2 c -fi=i e 

\/l0 \/l5 \/3 


2-P3/2 ! 2' 1 

2^3/2, |> | 

op 3 3 
Z ^l/2i 2> 2 

00 3 3 
^1/2, 2J 2 


a/To + 

x^5 e + 




3( 9 +2) V^(ff+2) n 

io c u 5^ Cu VTE 60 ° 

^(9+2) H(9+2) 2 (9-l) o 
5^ Cu 30 C u 3v ^ e U 

v^s-i) 2 (9-i) 9-4 
-fe 



(Table QUa) 



2-P3/2; !> I 



2P 



3 1 



3/2> 2' 2 



2P 



3 1 



1/2' 2> 2 



25 



3 1 
1/2' 2> 2 



2P 



1 1 



3/2' 2' 2 



2P 



1 1 



1/2' 2> 2 



25'l/2' |' I 



2P 



3/2 



2P 



3/2 



2Pi 



25 



1/2 
1/2 



5 3 
2' 2 



3 3 
2' 2 



3 3 
2' 2 



3 3 
' 2' 2 



2(9+2) 
5 e - 


9±1 P 
5 e - 


Vs 














lO-v/3 


11^(9+2) 
30\/3 


V8(g-i) 
3VT5 





v / io(9+2) 


V5(g-1) 
3V3 





x/30 


\/8(s-i) 

3^ 


V2{g-A) 





\/2(<?-i)„ 


3^ 














VE e - 








V3 e 



(Table flOb^ 



op 53 op, 33 
^3/2 > 2> 2 Z/ 3/2,2'2 



2 -Pl/2, |, | 



nc 3 3 
^l/2i 2' 2 




(Table flOh ) 



2-^3/2) 2' 2 



2 ^3/2, |, | 



2^1/2, f , ^ 



2^1/2) |) ^ 



2P: 



3/2 > 2' I 



2-Pl/2> 2 ' 5 



25i/ 2 , 2, | 



2^3/2) 1 5 ^ 
2^3/2, §, ^ 
2 A/2; 1 1 5 
2<S'i/2, §, 2 



2F 



3/2, I> 2 



2P n 



1/2' 2' 2 



25*1/2, |; 2 



9+2 „ 
9+2 „ 



.9-1 








eo 



9^ £ ° 


Vs(g+2) rii 

9 e ° 

10( g -l) , 



-eo 







_ 2 Cg-i) en 

o 

V2(g-4) 



9 



eo 



















en 





9 



18 e ° 
2y / 2(< ? -l) 



e 



9 






^eo 
^e 




e 



9 

9-4 „ 



e 









8«o 



(Table flOH) 



2-^3/2) 2> 



-\ 2P 3/2, |, 



•| 2P 1 /2,|,- 



2^1/2) |) 



25*1/2, 2' 



2^3/2, 2' I 
2 ^3/2, 2' I 
2-Pl/2) f ) ^ 
25*1/2 , |; 2 



2^3/2, §: 



9P 11 
Z/ l/2> 2' 2 

25*1/2, 2> 2 



3^(3+2) 
10 

^(g+2) 



10 
.9-1 



10 








10 

ll>/2(g+2) 
45 

4^(3-1) 
9\/5 





10(g+2l 
18 ' 







V10 

4^2(3-1) 
9\/5 

v^(g-4) r 
9 e 



v5(g-l) 
18 ' 

(3-4) - 
9 e - 









3 ' 








10(9+2) - 
18 

vg(g-ll p 
18 



5^(3+2) 
18 ' 

4(9-1) r 
9 e 







9 

(9-4) 
9 e 



4(9-1) - 
9 e 

V^(g-4) , 
18 1 











3 e - 






V2g 



e 



(Table HOb ) 



2-^3/2, f 5 



2Pq/9, I, 7? 2Pi/ 2 )|)| 25i/2,|,2 2P 3 / 2 ,7j,7j 2P 1 / 2 ,7j,?j 25x/2, 2J I 



3/2) 2' 2 



2-P3/2) 2' 
2-P3/2) |; 
2-Pl/2) f ; 
25*1/2, 2> 
2^3/2, 2> 
2^1/2, 2' 
25*1/2, 2' 



3^(3+2) 

— IS — e + 

_V2(3+2) 

10 e + 

>/l0 + 








^(9+2) „ 
10 e + 

11^(3+2) 

45 e + 
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Appendix B The adiabaticity condition 

In this appendix we discuss the conditions, which have to be satisfied in order to apply 
the general results of I to the n = 2 systems of hydrogen and deuterium. 

The condition to have a group of metastable and a group of fast decaying states was 
already discussed in sectional We found in (13.31) that we have to require \£\ < 250V/cm 
for \B\ = 0. Turning on the 3 field (\3 < 5mT|) one finds that the lifetimes of the n = 2 
states do not change substantially. 

Now we come to the question how slow the variation in time of the £ and B fields ought 
to be for the results of I to be applicable. In I we worked, for mathematical convenience, 
with a reduced time r and a total time T. We studied the limit T — > oo. 

Consider, for example, the two-state system of section IJH The leading terms for T — > oo 
are given in (I J4.25[) and (I J4.26I) and, with an estimate of the first correction terms, in 
(I JA.14I) and (I JA.15[) . The first correction terms should be small and thus we should require 

2c 12 



« 1 , (B.l) 



mm 



TAT 

4Cl2C21 «1, (B.2) 



T AT 

''21 



mm 

« 1 , (B.3) 
TAT v 1 

where ci 2 and c 2 i are defined in (I J4.20I) and (I JA.2j) . respectively. 

Now we revert to time t instead of the reduced time r. That is, we set T = To in the 
following, see (I J2 .3H . It is clear that the typical observation time will correspond to the 
lifetime of the 2S states, that is, we can set 

T « r s . (B.4) 

From the definitions (I J4.3H and (I J4.20H we then find 

C12 = T max |ai 2 (t)| , (B.5) 

0<t<T 

Mt)| = |(M|i^|2,t)| • (B.6) 

Now we insert for the state \l,t) the 2Si/ 2 , for the state \2,t) the 2P 1 / 2 states. We can then 
estimate the allowed speed of variation of the electric field using perturbation theory with 
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perturbing term D ■ 8£ in the mass matrix. We get in this way for constant magnetic field 

-—±—-—\2S 1/2 ,£(t),B) 
x (2S x ^I(t), B\ (-D) ■ (£(t + 8t) - Sit)) \2P 1/2 , S(t), B) , 



2P 1/2 , S(t + 6t), B) - |2P 1/2 , S(t), B) + 



%-\2P 1/2 ,£(t),B) = \\2S ll2 ,S{t),B)i2S 1 '^sXt),B\D 



dV 



V 



| «12 (*) | 



dSjt) 
dt 1 

1(25-!/^), B D • ^|2Pi/ 2 , 5 (t), B) 



er B (l) 



8S{t) 



Ot 



2P 1/2 ,S(t),B) , (B.8) 



(B.9) 



For hydrogen and deuterium we have (see (3.16) and (3.19) of Q]) 

er B (l) ^ 1 



Thus we get 



Vzs ' 



\a 12 (t)\ 



S = 477.3 V/cm. 



dS{t) 



VZS 



dt 



For the decay rates we have 



(B.10) 



(B.11) 



Ar 



min — 1 P 



S = 1 P 



Tp . 



[B.12) 



Inserting all this in (1B.5H and (IB.ip leads to the requirement 



2c 12 



TAr r 



1 1 

— max — 

Tp te[o,T] Sq 



Tp max — 

«6[o,n S 



dS{t) 



dt 
dS{t) 



dt 



< 1 



The same estimate is obtained from (IB. 31) . From (1B.2[) we get 



4ci 2 c 2 i 
TAr ■ 



TAr ■ 

J. Lil rnir 



2c 



12 



TAr r 



2c 



21 



TAr r 



rp max — 
te[o,T] £ 



95 (t) 



T max — 

te[o,T] S 



dS{t) 



Dt 



< 1 



(B.13) 



(B.14) 



Both conditions, ( IB. 13[) and ( IB. 141) . are satisfied if 



max — 

te[o,T] S 



dS{t) 



Ot 



Tp 



(B.15) 
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and 



max — 

te[o,n S 



dS(t) 



Of 



< 



T T S 



(B.16) 



Since we always have t$ ^> rp we see that ( IB. 161) already implies ( IB. 151) . 

Next we want to estimate the allowed rate of change for the magnetic field B(t). Similarly 
to (1B.9P we can estimate here 



|ai 2 (t)| 



±(2S 1/2 ,£(t), B(t)\n • ^Vi/ 2 , £(t), B(t)) 



(B.17) 



For the free 25* and 2P states - disregarding P violation - the matrix element in (1B.17|) 
vanishes. But we are considering states in an electric field where there is Stark mixing. For 
electric fields satisfying (13. 3p . \E\ < 250 V/cm, this mixing is at most 0(1). We shall thus 
estimate from (IB. 171) 



|ai 2 (t)| 



< 



L 



dB(t) 



Of 



In analogy to ( IB. 101) we define Bq by 



which gives 



o 



B n = 43.65 mT . 



(B.18) 



(B.19) 



^B.20) 



Repeating the same arguments as made above for the electric field we find the conditions 



max — 

te[o,T] B 



dB{t) 



dt 



and 



max — 

te[o,T] Bo 



dB(t) 



dt 



1 

T 



1 

rp 



(B.21) 



1 



(B.22) 



Here again ( IB. 221) implies ( IB. 211) since we have Ts 3> rp. 

We summarise the adiabaticity conditions for the electric and magnetic fields found as 
follows. The evolution of the — of course mixed — 2S states will decouple from that of the 
2P states if the electric field satisfies ( 13 .3p which guarantees for the lifetimes r$ ^> rp and if 
the rate of change of the electric and magnetic fields normalised to So and Bq, respectively, 
is smaller than the inverse lifetime r^T 1 of the 2S states. 
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Appendix C Geometric phases and flux densities 

In this appendix we will give the details of the calculations that lead to the results shown 
in sections |3] and HI In section IC.ll we consider a closed path in parameter space and 
transform the integral of the geometric phase into a surface-integral using differential form 
algebra. In section IC.2I we introduce geometric flux densities. In these sections we discuss 
general properties of geometric phases and, therefore, suppress the PV parameters in the 
labels of the states etc. In section IC.3I we use perturbation theory to identify the PC and 
PV contributions to the geometric phase and we derive the geometric flux densities that are 
used in section H] to visualise the geometric phase contributions. 



C.l. The geometric phase as surface-integral in parameter space 

Let R = (Ri, . . . , R r ) be the parameter vectors and C a closed curve in an r-dimensional 
parameter space 1Z. Let ^_{R) be a — in general non-hermitian — matrix function over TZ 
with non-degenerate eigenvalues E(a, R) (a = 1, . . . , N) for all R e C. 

The abelian geometric phase (I3.19P in parameter space reads (omitting the labels 8^2) 

laa{C) = j>dR-{^R\iV R \a,R) , (C.l) 
and can be written as an integral over a differential 1-form 

laa (C)=ij(^R\d\a,R) , (C.2) 
where the exterior derivative d is defined as 

i=l 

The exterior product (wedge product) of two 1-forms gives a 2-form. The wedge product A 
is antisymmetric, that is 

dRi A dRj = -dRj A dRi . (C.4) 

One can easily prove the following relations 

dRi A dRi = , (C.5) 
dd — , (Poincare lemma) . (C.6) 
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For an introduction to differential forms see 24| . As was already pointed out by Berry in his 



original work 



25], for a parameter space of dimension r > 3 one has to use the generalised 



Stokes' theorem to transform the integral ( 1C.2I) into a surface-integral. The generalised 
Stokes' theorem gives 

7 aQ (C) = i* (a^R\d\a, R) = i d(c^R\d\a, R) (C.7) 

JC=dF JT 

where T is an arbitrary surface in parameter space bounded by the curve C and lying fully in 
the regularity domain of d(a, R\d\a, R). From the Poincare lemma ( 1C.6I) and the asymmetry 
of the wedge product (1C.4j) we obtain for the integrand in (1C.7j) 



d(a, R\d\a, R) = (d(a, R\\d\a, R) 

= (d(c^R\) \(3, R) A (^R\d\a, R) 



- ((«» R W, R )) A (A R\d\a, R ) 



Here we have used the completeness of the eigenstates of ^f(R) 



(C.8) 



(C.9) 



and the following relations, which are easy to prove, 



(a, R\d\P, R) — — [d(a, R\j \0, R) , 

= ({aJt\d\a,R)) A (aJt\d\a,R) 



(CIO) 
(C.ll) 

Next we study the matrix elements (/3, -R|c?|a, R). For each parameter vector R we have 



(C.12) 
(C.13) 



R)\a,R) = E(a,R)\a,R) 
u7R\JL{R) = {o^R\E{a,R) 
(a = l,...,N) . 



Applying the exterior derivative to flC.12p . multiplying from the left with (/3,R\ and using 
the non- degeneracy of eigenvalues (which we supposed) we obtain for (5 ^ a 

{dM.R))pa 



((3, R\d\a,R) 



E(a, R) - E(p, R) 



(C.14) 
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with the 1-forms 

(fL#(-R))/3 Q = (dJOR)} \a, R) . (C.15) 

Using (E2D, (jOTj), (ELHD and flCH!) we finally get 



7UC) = if: [M^kAM*, (C = 3F;a = l,...,JV). (C.16) 



(E(a,R)-E(f3,R)) 



This is the geometric phase in the case of a non-hermit ian mass matrix, written as a surface- 
integral in parameter space. The result (1C.16j) is completely analogous to that in [25J] for 



hermitian Hamiltonians. However, in the case of a non-hermitian mass matrix, the geometric 
phase will in general be complex and therefore contribute to the exponential decay of the 
corresponding eigenstate. 

C.2. Geometric Flux Densities 

The mass matrix ^_(R) for an atom in an external electric (£) and magnetic (B) field 
has the form (see (T2T2B - ( 12^231) ) 

JL(E) = JU - D ■ £ -v- B , (C.17) 

where the parameter vector R is 

R=(R 1 ,...,R 6 ) = (S 1 ,£ 2 ,S 3 ,B 1 ,B 2 ,B 3 ) . (C.18) 

The exterior derivative of the mass matrix (1C.17j) is 



[R) = -D d£ - //• dB . (C.19) 

Thus, the matrix elements (1C.15P of d^£{R) are 

(d^_(R)) al3 = (<vR| (<L#)|/3, R) 

= -{c^R\D\p, R)-d£- (cvRH/3, R) ■ dB (C.20) 
= -D a p(R) ■ d£ - vjR) ■ dB . 
The 2-form occurring in (IC.16I) is 

3 

R)U A (dMR)h* = E (V-i^{R) d£i + ^ JR) dBi) 

U=i ' (C.21) 



A (D^R) d£j + {j, jJ3a (R) dBj) 
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Expanding the wedge product in (1C.21|) and using (1C.4I) leads to three different contributions, 



proportional to dSi A dSj, dBi A dBj and d£i A dBj, respectively. The geometric phase (1C.1 6j) 
can then be written in the form 



with the 2-forms 



(C) = [ ZS(H) + / 2ff(H) + / l£ B) (R) (C22) 

J T JT JT 



^J(R) = - E E Eh^BB^dSt^l iE, a ^ , (c.23) 

= 1 E E 2 - dft A dBL , (C.24) 

xm ^^-^^^ A ^ . (C25) 

For general external field configurations, the geometric phase (10.220 will get non-vanishing 
contributions from all three surface-integrals. 

Consider now a closed curve C in parameter space with constant magnetic field. If it is 
not prevented by singularities of the integrand we can choose also T to correspond to this 
constant magnetic field. Then only the integral over Xaa{R) will be different from zero. 
The 2-formXi« ) ( J R) (10231 is a product of two antisymmetric tensors. Introducing 

1 3 

dff ] = 2 E £ W dS > A dE i ' ( C - 26 ) 

i,;=i Wa {E{a,R) - E(P,R)) J 
(£ = 1,2,3) 

we can write (1C.23I) as a scalar product, 

^(R) = J£{R) ■ df £) = E J$ a (R)dfi e) , (C28) 

e=i 

and obtain from ( 1C.22I) 

7oa(C)= J J%(R)-dfW • (C.29) 

Obviously, <i/ ^ (10.260 is the oriented surface element in the 3-dimensional space of electric 
field strengths and J'^li.R) can be interpreted as a geometric flux density field in the space 
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of electric field strengths. The geometric phase for a given closed path C in the space of 
electric field strengths is given by the flux of the vector field ST^liR) through C. The 
graphical visualisation of the flux density field (1C.27I) turns out to be a very useful tool for 
choosing closed curves C in parameter space with desired properties, for instance minimising 
the decay rate for the metastable states along with maximising the geometric phase. 

For a closed curve C in parameter space with constant electric field the geometric phase 
laa{C) can be implemented in a similar way. We only have to replace in (1C.26I) to (1C.29|) 

(g\ (IS) 

£ by 13 , the 2- form Xad{R) by the 2- form X aa (R), and D. by /x. For general external field 
configurations where both £ and 6 vary we still can define the flux density fields J' ( fa{R) 
and J'^ai.R) an d write the geometric phase (1C.22|) as 

l aa (C) = I J<£{R) ■ df £) + I J<g{R) ■ dfW + / T^ a B \R) (C.30) 
Jt Jt Jt 

with lj^ 13 \R) from flC.25j) . However, visualisation in a single three dimensional parameter 

space is no longer possible. 



C.3. P-conserving and P-violating flux densities 

In this section we discuss the PC and PV contributions to the flux densities and geometric 
phases. Therefore, the dependence of quantities on the PV parameters is kept explicitly in 
this section. The mass matrix for the atom in external fields and including P-violation is 
given in (12.211) . For ease of notation we define 

5= {51 + 51) 1/2 , (C.31) 

2 c 

^p V = E i^pv . ( c - 32 ) 

8=1 

where 8^2 are given in (12.181) to (I2.20p . We can now write (12.211) and (1C.17I) in the form 

JLk$) = ^0 + ^pv > (C33) 
MR, 5)=^ + 5JZ_ VW -D £-v B (C.34) 

where ^ is the mass matrix for zero external fields and no PV contribution. We use 
perturbation theory to obtain a power series expansion in the small parameter 5 for the 
eigenstates of ^(R,5). The details of the formalism which applies to non-hermitian mass 

ri n 

matrices can be found in [7J, appendix C and [161 ]. appendix C. 
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The perturbation series expansion of the eigenvalues of ^(R, 5) reads 

E{a,R,5) = E(a {0 \R) + 0(5 2 ) . (C.35) 
There are no 5-linear terms in these eigenvalues due to time reversal symmetry invariance, 



see section 1.3 in 
of jM{R, S) reads 



161 ] . The perturbation series expansion for the left and right eigenvectors 



a, R, 5) = |a (0) , R) + R) + 0(5 2 ) , (C.36) 

^-PV,f3a(R) 

E(aW,R)-E(J3W,R 



\Q 



(1) .^EI^«) ^S iff . ^ 



a, i2, 5] = (a(°), R\ + 5(a«, R| + C(5 2 ) , (C.38) 

^PV,a/3(-R) 

R) — E(@(°\ R) 



[a 



In zeroth order we are left with the eigenstates \a^°\R) of the mass matrix j${R y 0) of 
f)C.34j) . In first order contributions proportional to the matrix elements of ^# PV occur 



which are defined as 

^pv, q/3 (^) = (^^Cr|^pvI/? (0) , R) • (G.40) 
The geometric flux density ( 1C.27I) which has PC and PV contributions, is given by 

(£=1,2,3) . 

The matrix elements of the dipole operator .D in (IC.41I) are 

D a/3 (R, 5) = (^R6\D\f3, R, 6) (C.42) 

and by using (10.36j) ff. we obtain, to first order in 5, 

D^iR, 5) = D™(R) + 6D%(R) , (C.43) 
D p a c JR) = {^Tr\D\^\ R) , (C.44) 



— a p {") - 2^ £7 ( a ( ) j _ E(j(°\R) 2-j E((3(°\ R) - E(<y(°),R 



(C.45) 
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Using flC.43j) to (1C.45I) . we are able to separate the PC and PV contributions to the 
geometric flux density (1C.41I) and find 



(£.pc). 



R) 



Dl%(R)D™ a (R) - (i «- j) 



{E{a(°\R)-E(fm,R)) 



(C.46) 
(C.47) 



-.PV 



PV 



iPC 



(E{aW,R) -E(p(°),R))' 



(C.48) 



The corresponding results for the flux density field ^T^a follow from straightforward 

substitution of £ by B and C by ^ in (fCTill to (jUHJ). 

We will now investigate the PC and PV character of the geometric flux densities ( IC.47j) 
and (1C.48|) . We define the quasiprojectors for the mass matrix ^_(R,0) = B, 0) 

fl034j) as 



Pi 0) (f,B) = \a w ,£,B)(a(°\£,B\ . 



(C.49) 



Here and in the following we write out R as (£,B), see (13.161) . With the completeness 
relation 



1 = 5>(°)(£,0) 



(C.50) 



we obtain from flC.471) and flC.441) 



p( 0) (£,B) 



xDjPg>(8,B)-(i~j) 



Consider now the parity transformation operator P defined by 



(C.51) 



P : x 



x 



(C.52) 



and its matrix representation V_ in the (n = 2) subspace of atomic states, 



V\2Lj,F,F 3 ) = (-l) L \2Lj,F,F 3 ) 



(C.53) 
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With V_ we have 

VxV} = -x , (C.54) 

P}V = 1 , (C.55) 

VDV} = —D , (C.56) 

V[J0 = v ■ (C.57) 

Thus, the mass matrix (10.340 satisfies for 5 = 

V_JL{£, B, 0) V} = V - D ■ £ - m • B) V} = JL{-E, B, 0) . (C.58) 

We will now use an argumentation analogous to the one in (7|, (3.43)ff. Assuming non- 
degenerate eigenvalues of the mass matrix for the considered field configuration, the quasipro- 
jectors ( 1C.49I) are the residues of the poles of the resolvent 

1 v R ( °\£,B) 



M£, B, o) - a E ( a(0) > £ ' B )-t 

where £ is an arbitrary complex parameter. Applying the parity transformation to the left 
hand side of flC39l we obtain with flC38l) 

p ! pt = 1 = V ^(-g'g) (C 60) 

Of course this must be equal to the parity transformation of the right hand side of (IC.59I) . 

and we get 

J- 2^ E f a (o) £ B) ^ E(a(°), -£, B) - £ ' { } 

Clearly, with our choice of B (see flA.lj) ) and our numbering scheme, table HI the states 
\a^°\£ = 0,B) are just the states \2Lj,F, F 3 ,£ = 0,B) which are eigenstates of V with 
eigenvalue (— 1) L . Turning on the electric field we find from continuity arguments that with 
the numbering scheme of table H] we have 

E(a {0 \ £, B) = E(a {0 \ B) , (C.62) 
VP£\£,B)V} = P£\-£,B) . (C.63) 

By inserting the identity (1C.55I) between each dipole operator and quasiprojector in (1C.51j) 
and by using flC36|) . §CM) and flC33|) we find 

j£/ c \£,B) = J% pc \-£,B). (C.64) 
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Thus, the PC geometric flux fields are indeed invariant under parity transformation of the 
external fields, 

P : (E,B) — > (-£,£) . (C.65) 

For the transformation properties of the PV geometric flux density fields we use 

V^p Y V} = -^ PV (C.66) 

and find with the same methods as for the PC flux density fields the result 

j£/ Y \£, B) = -Ji £ / V) (-£, B) . (C.67) 

Similar relations can be derived easily for the flux densities in magnetic field space, 

J^f G \S,B) = jff G \-S,B), (C.68) 
jff Y) {S, B) = -j£ py \-£, B) . (C.69) 

Finally, we consider the matrix elements of the time derivative 

^(R(t),5) = (aMf),S\^\(3,R(t),5) , (C.70) 

see (I3.20p . Inserting the expansions ( 1C.36I) and (1C.38I) we get 

where 

IkwiW)) = (a^R(t)\^°\R(t)) , (C.72) 



SwrtW)) = (aWMt)\^°\R(t)) + (a(^R(t)\^\R(t)) . (C.73) 



Now we consider a = (3. We have from ( 1C.37I) for all t 

{a^R{t)\a {l \R{t)) = (C.74) 

which implies 

(a(^(t)|^|a( 1 \i2(t)) + (a(^(t)|^|a«, J R(t)) = . (C.75) 



Inserting this in (1C.73j) for a = (5 and using (1C.37j) and (1C.39I) we get 

1 



E{a(°),R{t))-E{j(°\R{t)) 



d 



<^ PY ^(Rmi {0) ,R(t)\g- t \<y {0 \R(t)) 



(a^R(t) 1 1 17 (0) , R(t) )l PV , 7Q (fl(t) ) 



From 



[a 



W,R(t)\^°\R(t))=5, 



cr; 



we get 



(a^R(t)\^°\ R(t)) + (a(^(i)|4l7 (0) , «(*)) = . 



'at 1 



Inserting this in ( 1C.76j) and using ( 1C.72I) leads to 

1 



&PV,a*(R(t)) = E 



£(a(°),tf(t))-£(7(°),H(t)) 
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(C.76) 



(C.77) 



(C.78) 



(C.79) 



X <^ ^P V , a7 WMp C , 7 a(^(0) + ^PC, Q7 (#(0)^PV, 7 aW)) 



Splitting up 5@py aa (R(t)) into the contributions proportional to Si and 5 2 we write 

«&v« W)) = Si@£La(R(t)) + • (C80) 



Inserting here fl079l leads to the results (I3T2D and fl3T25|) . 
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